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SUMMARY 


SYMBOLS 


The theory of Lagrangian muUipliera is applied to the prob- 
lem of finding both upper and lower limits to the true compressive 
buckling stress of a clamped rectangular plate. The upper and 
lower limits thus bracket the true stress, which cannot be exactly 
found by the differential-equation approach. The procedure for 
obtaining the upper limit, which is believed to be new, presents 
certain advantages over the classical Rayleigh-Ritz method of 
finding upper limits. The theory of the lower-limit procedure has 
been given by Trefftz but, in the present application, the method 
differs from that of Trefftz in a way that makes it inherently more 
quickly convergent. It is expected that in other buckling prob- 
lems and in some vibration problems the Lagrangian multiplier 
method of finding upper and lower hmiis may be advantageously 
applied to the calculation of buckling stresses and natural 
frequencies. 

INTRODUCTION 

Many important problems that cannot be exactly solved by 
the differential-equation approach and must therefore be 
analyzed by approximate methods arise in the buckling and 
vibrations of thin plates. The theory of Lagrangian multi- 
phers can be a powerful tool in the analysis of many of these 
problems. The present paper presents the details of applica- 
tion as well as the fundamental principles of the Lagrangian 
multipher method by demonstrating the use of the method 
to obtain both upper and lower limits to the true compressive 
buckling stress of a rectangular plate clamped along all edges. 

The procedure for obtaining the lower limit is similar to a 
method used by Trefftz (reference 1) and recently described 
by Reissner (reference 2). The present lower-limit method 
differs from that of Trefftz, however, in a way that makes it 
inherently more quickly convergent. The upper-limit pro- 
cedure, which does not appear to have been presented pre- 
viously, is simpler than the usual Rayleigh-Ritz method and 
may be expected to permit the computation of more accurate 
results with less labor. 

In a recent treatment of the problem of compressive buck- 
ling of clamped plates, extensive calculations of lower limits 
were made by Levy (reference 3) by means of a procedme 
equivalent to the Trefftz method. The results were esti- 
mated by Levy to be within 0.1 percent of the true results. 
In order to Dlustrate the methods of the present paper, upper 
and lower limits to the buckling stress of a square plate are 
computed to within 0.1 percent of each other; a positive 
check on the accuracy of Levy's results is thus obtained. 


a length of plate, in direction of stress 

b width of plate, peipendicular to stress 

0 aspect ratio (a/b) 

t thickness 

ft Poisson's ratio 

E Young's modulus of elasticity 


D plate stiffness in bending 


X plate coordinate in direction of stress 

y plate coordinate, perpendicular to direction of 

stress 

w plate buckling deformation, normal to plane of 

the plate 

(T* critical compressive stress, in z-direction 

k critical compressive stress coefiBcient in the 

formula o-,=Ar 

V internal energy of deformation 

T external work of applied stress 


a„ Fourier coefficient of cos — - 

" a 

br Fourier coefficient of cos 

a 


a^n Fourier coefficient of cos — - cos — 

a 0 

i, j,m,n,p,q even integers 
r, s odd integers 

Kronecker delta (1 if m=n; 0 if m?^n) 

A I 

Of, X, hj, fit, 7 Lagrangian multipliers 

THEORETICAL BACKGROUND 

Eayleigh-Eitz method. — The Rayleigh-Ritz energy method 
for determining the critical stress of a thin plate consists of 
the following steps: 

(1) The deflection surface of the buckled plate is expressed 
in expanded form as the sum of an infinite set of functions 
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having undetermined coefficients. In general, each term of 
the expansion must satisfy the geometrical boundary condi- 
tions of the problem. 

(2) The energy of the load-plate system is computed for 
this deflection surface and is then minimized with respect to 
the undetermined coefficients. 

(3) This minimizing procedure leads to a set of linear 
homogeneous equations in the undetermined coefficients. 
These equations have nonvanishing solutions only if the 
determinant of their coefficients vanishes. The vanishing of 
this stability determinant provides the equation that may be 
solved for the buckling stress. 

When the set of functions used is a. complete set capable of 
representing the deflection, slope, and curvature of any possi- 
ble plate deformation, the solution obtained is, in principle, 
exact. Since, however, the exact stability deteiminant is 
usually infinite, a finite determinant yielding approximate 
results is used instead. 

lagrangian multiplier method. — ^The Lagrangian mul- 
tiplier method follows the general procedure outlined for 
the Rayleigh-Ritz method, with but one outstanding change. 
The restriction in step (1) that the boundary conditions be 
satisfied by every term of the expansion is discarded and is 
replaced by the condition that the expansion as a whole sat- 
isfy the boundary conditions. This condition is mathemat- 
ically satisfied in step (2), during the minimization process, 
by the use of Lagrangian multipliers. 

The fundamental advantage of the Lagrangian multiplier 
method lies in the fact that, with the rejection of the neces- 
sity of the fulfillment of boundary conditions term by term, 
the choice of an expansion is much less restricted. In the 
clamped-plate compression problem, a simple Fourier expan- 
sion may be used instead of the compheated functions 
assumed in the Rayleigh-Ritz analyses of this problem (refer- 
ences 4 and 5). Furthermore, the orthogonality properties 
of the simple Fourier expansion lead to energy expressions 
of a simplicity that is instnunental in permitting accurate 
computations. 

Approximate solutions: upper and lower limits. — 

The Lagrangian multiplier method, as well as the Rayleigh- 
Ritz method, gives a theoretically exact solution for the 
buckling stress; but ordinarily only approximate results are 
obtained because, of tlie practical necessity of considering 
finite rather than infinite determinants. In the Rayleigh- 
Ritz method the approximate result is always higher than 
the true buckling stress. In the Lagrangian multiplier 
method, however, it is possible to obtain approximate solu- 
tions in two different ways, wliich permit the computation 
of a lower limit as well as an upper limit to the time buckling 
stress. As determinants of higher order are used to obtain 
approximations of hi^er order, both the upper-limit and 
lower-hmit results approach the true buckling stress. Thus, 
the Lagran^an multiplier method can be used to provide a 
result to within any specified degree of accuracy. It may 
be expected, furthermore, that a particular determinant in 
the Lagrangian multiplier method ought to yield a more 
accurate result than a determinant of equal order in the 
Rayleigh-Ritz method. 


LAGRANGIAN MULTIPLIERS 

The procedure used in applying the fundamental mathe- 
matical principles of Lagrangian multipliers is described in 
this section; a general proof of the validity of the method is 
given in the appendix. 

Let it be required to minimize a function of N variables 

/{XuXt,X3,...Xif) ( 1 ) 

where the x’s are not independent but are bound together 
by the relationship 

p{Xi, Xi, Xi, . . . Xff )=0 (2) 

Lagrange’s method of simultaneously minimizing / and 
satisfying the constraining relationship (2) is to minimize 
the fimction 




with respect to the x’s. The quantity X is the undetermined 
Lagrangian multiplier. The necessary conditions for min- 
imizing/ then become 


(K=,,2,3,...A0 

p=0 (2) 

Note that these expressions are iV-t- 1 equations in the N -[- 1 
unknowns Xi, Xt, ... Xs and X. 

If there are two relationships that constrain the «’s; 
that is, if 

. . ■ a:w)=0 


— 0 


two Lagrangian multipliers are then needed. The function 
to be minimized becomes 


/ — 


and the minimizing equations are 



dxs: 


dzfc 


(K= 1 , 2 , 3 , ... AO 


Pi=0 


<Ps=0 


The method is easily extended to cover the case of any 
number of constraining relationships. 

PRELIMINARY ILLUSTRATIVE EXAMPLE 

Before the main example is given, a simpler buckling 
problem will be analyzed by the Lagrangian multiplier 
method in order that the method of application of Lagrangian 
multipliers may be most clearly presented without the ob- 
scuring details of analysis of more complicated problems. 
This elementaiy problem requires tlie use of but a single 
Lagrangian multiplier, which leads to a single stability 
equation. 
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Consider a square plate, clamped along two opposite 
edges, simply supported along the other two edges, which is 
loaded in compression on the simply supported edges. (See 
sketch.) From the exact solution of this problem, the deflec- 
tion surface of the plate is known to be sinusoidal in the 
a:-direction. The deflection in the i/-direction, known to be 
symmetrical, must satisfy the clamped-edge boundary 
conditions; that is, zero deflection 

w(x,0)=w{x,a)=0 (3) 

and zero slope 

|?(x,0)=|^(a:,a)=0 (4) 

The present method uses a cosine-series expansion, whereas 
the Trefftz procedure would use a sine-series expansion. The 
problem is solved by both methods for comparison. 



Cosine-series solution. — In the cosine-series solution the 
expansion 

m 

. mvx 71/icy .... 

«j=sm— > a„cos-^ (5) 

may represent the deflection surface having m half-waves 

in the x-direction, since the Fourier series of even cosines is 
a complete sy mm etrical set. 

The boundary conditions (equation (4)) on the slope are 
satisfied by each term of the expansion; however, in order 
that w satisfy the conditions of equation (3) on the edge 
deflection, it is necessary that 

00 

^ a,=0 (6) 

Equation (6) is a constraining relationship on the a's and 
as such win be introduced in Lagrange's minimization process. 

As in the Rayleigh-Ritz method, the internal energy of 
deformation and the external work of the stresses are then 
calculated. Using the value for w as given by equation (5) 


in the general foimulas (reference 6, equations (199) and 
(201) (modified)) 


gives 


(7) 


( 8 ) 


V (m*-i-7iT(l + Oa, 


Zj 

n=0,2,4 

^ (l + 5on)a»* 

The usual Rayleigh-Ritz procedure requires that the 
expression 

V—T (9) 


be a minim um with respect to the a’s. In the present ex- 
ample, however, the a’s are not independent but are bound 
by equation (6). Hence, mathematically stated, the ex- 
pression V — T must be a minimum, subject to the constraint 
relationship on the a’s . _ 

CO 

^ a.=0 (6) 

»ioi2,4 

Solving this minimiz ation problem by Lagrange’s method 
makes it necessary to minimize 

rr 

(H-D-X^a» (10) 


with respect to the a’s. 'The necessary conditions for a 
minimum then become 


d(V-T) 

dOf 



(i=0,2,4,...) 





a»=0 


or, upon differentiation and simplification, 

4/7® 

(l+U[(m*-f-j*)*-m*t] o,-^X=0 (i=0, 2, 4. . . .) 0 1 


2 «»=0 ( 0 ) 
n-OA^ 
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Solving equation (11) for % and substituting into equa- 
tion (6) gives tbe stability equation that detennines k: 


1-0,2, 4 

For a particular number of half waves m, this equation 
may be solved by evaluating the series for several trial 
values of k and interpolating to find the k that makes the 
series vanish. The correct value of m is that which gives 
the lowest value of k. For two half waves (m=2) in the 
loaded direction, the theoretically exact value of A=7.69 
(reference .6, p. 345) is obtained when only 10 terms of 
equation (12) are computed. 

Sine-series solution (Trefiftz method). — ^The same problem 
will be treated in the manner suggested by Reissner (refer- 
ence 2), which is siniilar to Trefftz's method (reference 1). 

Let 

fl» 

. wixa; V' , . rvy 

w=sai — (13) 

T '^,6 

The boundary conditions on deflection (equation (3)) are 
now satisfied term by term, but the conditions on the edge 
slopes (equation (4)) are satisfied only by making 



(14) 


Now, the expression V— T is computed from formulas (7) 
and (8) by using the value of w given by equation (13) j 
then, by application of Lagrange’s procedure, 


00 

{V-T)-y^ rh 


(15) 




must be minimized with respect to h, where s=l, 3, 5, . . .. 
The Lagrangian multiplier is y. 

The minimization equations are 


stability equation, the Lagrangian multiplier method is 
preferably used to satisfy the zero-deflection condition rather 
than the zero-slope condition. Slope is the derivative of 
deflection, and, in general, differentiation of a Fourier scries 
makes it more slowly convergent. 

THE COMPRESSIVE BUCKLING OF A RECTANGULAR PLATE 
CLAMPED ALONG ALL EDGES 

The previous elementary example required only a simple 
Fourier expansion and but-xine Lagrangian multiplier to sat- 
isfy the boundary conditions. The more difficult problem 
of finding the buckling stress of the rectangular plato 
clamped on all edges and loaded as shown in the accomi>any- 
ing sketch necessitates a double Fourier series, as well as an 
infinite set of Lagrangian multipliers to satisfy the boundary 
conditions. 



Boundary conditions. — The boundary conditions of the 
problem are: 

Zero deflection, loaded edges 

w(0, y)=v){a, y)=0 (18) 




4a=' 


7«=0 (5=1, 3, 5, . . .) (16) 


2 rh=Q (14) 

Solving equation (16) for b, and substituting in equa- 
tion (14) gives as the stabihty equation 


CO 

E 

>- 1 , 8,8 




(17) 


Comparison and discussion of results. — The series 
in equation (17) converges approximately as 1/s*, whereas 
the series in equation (12) converges approximately as 1/j*. 
Because of the more rapid conver gen ce obtained in ^e 


Zero deflection, unloaded edges 

w{x, 0)=iy(ar, b)=0 (19) 


Zero slope, loaded edges 

^(0.s)=g(«,y)-0 (20) 

Zero slope, unloaded edges 

^{x,0)=^{x,b) = 0 ( 21 ) 

Fourier expansions. — In order to achieve a rapidly 
convergent solution, the principles established by the pre- 
ceding e^mple are used as the basis for choosing the Fourier 



LAGRANGIAN MULTIPLIER METHOD OF FINDING UPPER AND LOTSTER LIMITS TO CRITICAL STRESSES 


217 


expansioDs to satisfy, term by term, the conditions of zero 
slope rather than those of zero deflection. 

The buckling deformation corresponding to the lowest 
buckling stress is always symmetrical perpendicular to the 
direction of load but, depending on the aspect ratio of the 
plate, may be symmetrical or antisymmetrical in the direc- 
tion of load. Thus, for symmetrical buckling, let 


Since this Fourier series must vanish, each infinite series that 
constitutes a coefficient of a cosine term must vanish. (All 
the Fourier coefficients of the Fourier expansion of the func- 
tion zero are zero.) Hence 

CO 

^ a„j=0 C/=0,2,4, . . .) (25) 

m=Q,2,i 


at m 

E \\ mvx nvy 

> cos -^cos-^ (22) 

m-0r2,4 n=^4 

and, for antisymmetrical buckling, let 

^ 2 ^ ^ 

It is sufficient, for purposes of demonstration, to consider 
only the case of sjunmetrical buckling. Hereinafter, w there- 
fore refers to the value given by equation (22). 

Energy expressions. — Using the expansion given by 
equation (22) in the evaluation of the general energy and 
work integrals of equations (7) and (8) gives 

at at 

m-0,2,4 fl-0,2,1 


By expressing the fact that there is zero deflection along 
the unloaded edges (equation (19)), it can be similarly shown 
that 

eo 

^ (f=0,2,4, . . .) (26) 

tt^or2,4 

Now, Tmustbe a minimum with respect to the a% 
which are boimd by equations (25) and (26). As the prob- 
lem now stands, however, it is not in the form to which 
Lagrange’s minimization process can be applied since V— T_ 
does not contain Ooo, whereas the constraint relationships 
do contain Ooo. Hence, is eliminated from the con- 
straint relationships by subtracting the first of equations 
(26), the equation for f=0, from the first of equations (25), 
the equation for j—Q. The final set of necessary con- 
straining relationships on the minimization of the energy 
expression (24) 


T 


CO CD 

8a Zj Zj 

m=0, 2, 4 saO, 2,4 




Then 

V-T=^ 

where 

-^=[(mH7iW^-*m*^’](l-|-5,»o-l-5o.). 
Note that V— T is independent of a«ij since ^=0- 


2 2) 


(24) 


m-0,2,4 n-0,2,4 


Constraining relationships. — The boundary conditions 
of zero slope (equations (20) and (21)) are satisfied by 
each term of the expansion of equation (22), but the condi- 
tions on deflection (equations (18) and (19)) must be satisfied 
by the expansion as a whole. Substituting winto equation (18) 
gives, along the loaded edges. 


at CD 

W(0, y)=V){a, y)==^ cErt)+cos ^ ^ 


ffl.-0,2,4 


ia-0,2,4 

4xy 

jb’^ 2,4 


+cos^^ a„ri 


-hcos^^ a,^ 

+ . 


in-0,2,4 

=0 



then becomes 


CO 

flmO — Oo,= 0 

09 

a«/=0 

;^0r2,4 
00 


m-2, 4,6 n=2,4,6 


(i=2,4,6,...) 


y (27) 


0 (1=2,4, 6,...) 


n-0,2,4 


Theory of upper and lower limit solutions. — ^A theoretically 
exact solution to the problem would be obtained if the 
energy expression (equation (24)) were minimized with 
respect to all the a’s and at the same time all the rela- 
tionships (27) were satisfied. This result follows from the 
facts that: (a) the expansion of equation (22) is a complete 
symmetrical set, capable of representing the exact symmet- 
rical buckling deformation, and (b) the fulfillment of the 
conditions of equations (27) ensures that the boundary con- 
ditions are completely satisfied. An exact solution is not 
possible, however, because it would involve an infinite deter- 
minant, so that two different modifications of the ideal pro- 
cedme are used to obtain approximate results. One of these 
methods gives an upper limit to the true buckling stress, 
whereas the other gives a lower limit. 
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An upper limit to the buckling stress can be found by 
arbitrarily setting some d's equal to zero, minimizing expres- 
sion (24) with respect to the remaining a’s, and satisfying all 
the constraint relationships (27) . An upper limit is obtained 
inasmuch as arbitrarily setting somff~of the Fourier co- 
efficients equal to zero has the effect of restraining the deflec- 
tion of the plate, which in effect stiffens the interior of the 
plate and increases the stress required to buckle it. 

A lower hmit to the buckling stress can be found bj’- 
minimizing expression (24) with respect to all the a’s but 
satisfying only some of the constraining relationships (27). 
Neglecting some of the constraining relationships has the 
effect of giving tlie plate greater freedom at the edges and 
hence reducing the stress required to buckle the plate. 

Lower limit solution. — In accordance with the require- 
ments for a lower limit, the constraining relationships (27) 
will be satisfied only up to j=q and i=p. By Lagrmige’s 
minimization process, the function to be minimized is then 


(?= 


ir^Db 


O 03 / 09 CO \ 

j-O.2,4 n-0,2,4 \fn-2,4,6 n-2,4,6 / 


q 09 p 09 


^-2,4,6 m=0,2,4 


i-2,4,6 n-0,2,4 


The a, Hi . . . /i, are Lagrangian multipliers. 

The equations for minimizing V— T with the constraining 
relationships (27) on the a’s then become 


bG 

Sflins 


0 


(m,n=0, 2, 4, 



Equations (27) up to i=p,j=q 


(28) 


Evaluation of bGjbamn gives 

3^ ®«n~ «(5o»— 5mo) — X»— /Xm=0 (29) 

where X„ appears only if 2^n^q and jum appears only if 
2^m^p. From equation (29) 


< j+e — 0 (d, e — 2,4,6, . . .) 

Hence for any particular d, e, either 


or 


0 


«+« 


The first alternative, however, ordinarily would require k to 
be very high, corresponding to the budding stress of a 
buckling mode with many waves in both directions. For the 
low'est buckling load, then. 


flMit. «+«=0 (d, e=2, 4,0, . . .) 

It is therefore necessary to be concerned with only the other 
o’s, which, from equation (29), are 


_ 4a^ j /X _i_ \ 

4cd^ 


(m,n?^o) (30) 


In equations (30), does not appear if ft>gr and does 
not appear if m>p. 


Substituting the values of a’s given by equations (30) back into the constraining relationships (27) up to j=q, i=P gives 

P 


\m-2,4,6 n-2,4.9 / m-2,4,8 n-2,4,6 

— Anja+ ^ ^ Am}=0 (j=2,4,6, . . . g) 


m-2, 4, 6 


m =0.2.4 


00 (j 

A<bXb= 0 (i=2^4, 6, . . . p) 

n-l),2,4 n=2,4.6 


(31) 


These equations form a set of i(p+g) + l linear homo- 
geneous equations in a, /i. ... X* ... X,. Since when 
buckling occurs the a’s are not all zero, by equations (30), 
the Lagrangian multipliers are not all zero. In order that 
equations (31) be compatible, the deteraiinant of the coeffi- 
cients of the Lagrangian multipliers must vanish. The 
vanishing of this stability determinant provides the deter- 
minantal equation that may be solved for k by substitution 
of trial values and interpolation. 


That certain elements of the determinant consist of an 
infinite series of terms is evident; these series converge 
rapidly. Since such rapidly convergent series are cal- 
culable to-any degree of accuracj', they may be considered 
as known quantities. Each value of represents the 
potential-energy contiibution of a term in the expansion for 
w; hence, the effects of infinite subsets of expansion terms 
enter into this solution. Thus, for p=g=2, the expansion 
terms corresponding to the a’s showm in table I enter into 
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the solution; similarly, the terms represented in table II 
enter into the solution ^rhen f=q=A.. 


TABLE I.— FOUR INFINITE STRIPS OF FOURIER 
COEFFICIENTS OF EXPANSION TERMS 


w 

n 

0 

2 

4 

e 

8 

. . . 

0 

0c« 


att 

au 

a» 


2 

On 

an 

an 

an 

tta 

... 

4 

Oti 

On 





6 

a» 

<in 





8 

Ott 

au 






TABLE II.— SIX INFINITE STRIPS OF FOURIER 
COEFFICIENTS OF EXPANSION TERMS 


IB 

» 

2 

4 

6 

8 

• ■ - 

0 

ac« 


041 

an 

osi 


2 

on 

aa 

041 

an 

ou 

. . . 

4 

004 

au 

044 

an 

au 

. . 

6 

an 

aa 

041 




s 

an 

an 

040 





Upper-limit solution. — The lower-limit solution satisfied 
only some of the constraining relationships (27) but assumed 
the existence of all the Fourier coefficients. If an upper 
limit is to be obtained, it will be necessary to satisfy aU the 
constraining relationships while arbitrarily assuming some 
a’s to be zero. 

As a direct result of the necessity of satisfying all the 
constraining relationships in the upper-limit solution, it is 
found that the first of equations (27) is redundant and may 
be discarded, since it is automatically satisfied when aU the 
remaining equations (27) are satisfied. As a proof of this 
redundancy, the conditions 

«» 

. (j'=2,4, 6, . . .) 

m-0,2,4 

are summed over j and subtracted from the sum of the 
conditions 

e» 

^ a*,=0 (i=2,4,6, . , .) 

■ 

over i and give 



It is to be emphasized that the redundancy of a constrain- 
ing relationship is a peculiarity of only the upper-limit solu- 
tion, since, as shown by the proof given, the redundancy 
depends on the fact that all the constraining relationships 
must be satisfied. 

With the elimination of the redundant condition, the 
necessary constraint relationships become 

2 0=2, 4, 6,...) (32) 

nt^0j^4 

e» 

^ a,„=0 (i=2,4,6,...) (33) 

n-oTl4 

At this point, in accordance with upper-limit theory, it 
is necessary arbitrarily to set certain a's equal to zero. It is 
possible to take advantage of the Lagrangian multiplier 
method by allowing infinite rather than finite sets of a's 
to exist and still to obtain a stability determinant of finite 
order. Thus, infinite strips of coefficients of the type shown 
in tables 1 and 2 can enter into the solution. In the lower- 
limit case, the existence of aU coefficients was assumed, but 
the coefficients Up+i, s+, were proved to be zero; in this 
upper-limit solution, it will be arbitrarily assumed that 
these same a’s are zero; thus 

8+«=fi (c^) c=2, 4, 6, ... ) 

The constraining relationships (32) and (33) become 


Clmi — 0 

(i=2,4, 6, ... 2) 

(32a) 

a„)=0 

(J=3+2, 2+4, . . .00) 

(32b) 

at„=0 

(i=2, 4, 6, ... p) 

(33a) 

atn=0 

(i=p+2,p+4, ... 00 ) 

(33b) 


n=0, 2, i 

The function to be minimized is 


ni=0,2,4 a=0,2,4 1-2, 4, 6 m-0,2,4 ■ 

to p p to 

“ E E E E ® 

1-9+2, e+4 tb- 0.2,4 J-2,4, 6 n=0, 2, 4 

- V 

2j 2j 

I-P+2.P+4 n=0, 2, 4 


1-9+2, 9+4 7B-0.2,4 4-2,4, 6 n=0, 2, 4 
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The first double summation of this equation extends over 
only the values of m and n such that 

if n^q_ 


Setting 
to exist, 


71^2 if 

0 then gives, for aU the a's arbitrarily allowed 


— 




/^m)] 


where Xo and hq do not exist. Substituting back into the 
constraint equations (32a), (32b), (33a), and (33b) gives 


flO 

Amj+ 

00 

:0 (i=2,4,6, . . 

• 2) 

(34a) 

m=0r 2, 4 


m-2,4,6 




P 

Amy+ 


:0 (i=2+2,2+4,. 


(34b) 

m-0,2,4 


ni-%4,6 




E 

A-iaXa+Alt 

^ Atn=0 

(i=2, 4, 6, . . . 

P) 

(35a) 

n-2^4,6 


»=0^2,4 




2 

.^i«Xn+ yt 

^ A,.=0 

(t=p+2,p+4, . . 

.co) 

(35b) 

n-2,4,9 


B=0, 2 , 4 





These equations involve all the Lagrangian multipliers. 
They can be reduced to a set of equations, however, in 
tii ... m. the following manner: 

From equation (35b), for i=p+2, p+4, . . . 




y-i- 






n-0,2,4 


From equation (34b), for j=2+2, 2+4, • . . 




n;-0i%4 


Substituting these expressions for ni and X/ into equations 
(34a) and (35a), respectively, gives as the final stability 
equations: 




m>-0,2,4 m»2,4,6 



=0 (i=2,4,6, ... 2) 


q m 


I y ^ I 


■^mn 


(30) 


eo _ 

^in=0 (i=2, 4, 6, 


P) 


n-0,2,4 


Equations (36) form a set of ^ (p+2) linear homogeneous 

equations in X2 . . . Xg, /t2 . . . /*,. The stability determinant 
is the determinant of the coefficients of the X's and ft’s. 

It is of interest to note that in the usual Rayleigh-Ritz 
solutions only finite sets of expansion terms are ever taken 
into account, and the order of the determinant obtained is 
ordinarily equal to the number of terms considered. It is 
then reasonable that a particular determinant obtained by 
the Lagrangian multiplier method, which considers infinitely 
more expansion terms than a Rayleigh-Ritz determinant of 
equal order, may be expected to give a more accurate result. 

Numerical example. — ^For the case of a square plate, 18= 1, 
upper and lower limits were computed. The results for 
the buckling-stress coefficient k were:. . 


Approximation 

Lower limit 

Upper limit 


9.99 

lau 

ia.08 

Second; p=g=4. 

- 

— 

10.07 


The expectation that the Lagrangian multiplier method 
should give closer upper limits than the Rayleigh-Ritz 
method, for a given-order determinant, can be confirmed for 
this example. A second-order Lagrangian multiplier deter- 
minant gives an upper limit of k= 10.11, whereas Maulbetscb 
(reference 4) and Smith (reference 5) use complicated de- 
flection functions in the Rayleigh-Ritz method to derive 
third-order determinants that give, respectively, 4=10.48 
and 4=10.11. 

It is seen that the second approximation, requiring the 
evaluation of a fourtli-order determinant for the upper limit 
and a fifth-order determinant for the lower limit, definitely 
establishes the value of 4 to within 0.1 percent. 

Levy (reference 3) used an ingenious method of obtaining 
lower limits that is, in fact, equivalent to the Trofltz method 
of using double sine series and satisfying the zero edge-slope 
condition by the Lagrangian multiplier method. On the 
basis of computations involving determinants up to order 
twenty. Levy concluded that his results obtained from tenth- 
order determinants are within 0.1 percent of the true results. 
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inasmuch as Levy obtained ^:= 10.074 for the square plate, 
the present relatively simple upper- and lower-limit calcula- 
tions show that his estimated limi t of error is correct for 
this case. 

CONCLUSIONS 

1. The Lagrangian multiplier method can be used to 
compute accurate upper and lower limits to the compressive 
bucklipg stress of a clamped rectangular plate, thereby 
bracketing the true result. 

2. From a consideration of rapidity of convergence toward 
the exact solution in clamped-plate problems, it is preferable 
to use an expansion that satisfies the zero-slope boundary 
conditions term by term rather than the zero-defiection 
boundary conditions. 

3. Because of the fact that the Lagrangian multiplier 
method permits the effects of infinite subsets of expansion 


terms to enter into the solution, it is believed that a particular 
stability determinant derived by the Lagrangian multiplier 
method will, in general, yield a closer upper limit than that 
obtained from a determinant of equal order in the Rayleigh- 
Ritz method. 

4. It is expected that the method of Lagrangian multipliers 
may be useful in the analysis of other stability and vibration 
problems. In particular, the method may be immediately 
applied to the determination of vibration frequencies of 
clamped plates, and to the determination of buckling stresses 
of clamped plates under compression in two directions. 


Langley Memorial Aeronautical Laboratory, 
National Advisory Committee -for Aeron.4.utics, 
Langley Field, Va., May S, 1946. 



APPENDIX 

GENERAL PROOF OF THE METHOD OF LAGRANGIAN MULTIPLIERS 


Let it be required to minimize 

f{xi, xj, Xi, . . . Xjv) (Al) 

where the N I's are bound by the P independent relation- 
ships (P<iV) 

<Pj (2i) 3%) • • • 3;^r) = 0 (t7=l, 2, 3, ... P) (A2) 

It will be proved that the equations for determining the 


minimizing values of the x’s are: 


lx lx Sv’3 


■ ■+"'5S-“ (X=1.2,3,...A0 (A3) 


Xs, Xi,... Xjv)=0 (J=l, 2, 3, ... P)' 


(A2) 


The X’s are Lagrangian multipliers; these N+P equations 
deteimine N x’s and P X’s. 


If the values of only N—P x’s are known, the remaining P x’s are determined from the P relationships (A2). For 
convenience, consider the last P x’s in equation (Al) to be dependent upon the first N—P x’s. Then for / to be a 
minimum its fihst partial derivatives with respect to the independent x’s must vanish, or: 


a/ . a/ axjT-p+i , a/ 

axjf ax^_p+i axjif ax^_p+a ax_jf 


(M=l,2,3, . . . (A"-P)) 


(A4) 


But each of these equations contains P quantities that cannot be directly evaluated — the derivatives of the depend- 
ent variables with respect to the independent variables. For each value of il/, these P derivatives are determined by 
differentiating each of the P constraint relationships (A2) with respect to x^. Thus, 


bifij 


b<pj axjv-p+i 1 axjv_f >+2 1 I 




=0 


(J=l,2,3, . . . P) 


(A6) 


Now, for each particular value of M, equation (A4) and the P equations (A6) make up a set of P+1 linear homo- 
geneous equations in the P+1 quantities 1, • 


^ gjj^gg jjjggg quantities are surely not all zero, 


■<'3ir 


dx 


■Af 


the determinant of their coefficients must vanish. Hence, it is found that for / to be a minimum it must 
be true that: 



ay 

. 



£>xjf 


^X;vr-P+s 

dxjif 


dv5i 

b<pl 

i><pi 



&Xif 


dx^r_p+2 

dXjjr 


dv?a 

b<pi 


d<p2 

=0 (A/=l,2,3, 

iiXM 

^X^_P+l 

dXjff_P4.2 

dXjir 




3<£>p 

d^P 


3x^f 


dXy_p+2 

2>xjv 



necessarily 


(A6) 


It will now be demonstrated that these necessary minimization equations will hold if eauations (A3) liold. Interchanging 
the rows and columns of the determinant in equation (A6) gives: 


ay 


b(p2 

^<Pp 

/O 

dxjif 

^Xif 

bxsf 

^y 

d(Oi 

b<P2 

^<pp 

dx^_P+i 

dx^r_p+i 

^Xjv_p+i 

dXAT_i>+l 

¥ 

b<pi 

5^02 


dXAT_p+2 

Sx^f _P^ 2 

dxj\r_p.(.2 

dX^^_P4.2 

¥ 

d^i 


^<Pp 


bxtf 

bXif 

dxw 


(.1/=1,2, 3, ...(A^-P)) 


(A7) 
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The vanishing of this determinant is, however, precisely the condition of compatibility of the equations 




+ +..■+ -0 


- ¥ 
SaTiv-p+i 


^<Pi 

* ^w-p+i 




+Xi . . ■ +Xp5-!^=0 


btfip _ ^ 
&3:jif-p ^-1 


5?^+>'-5r®-+>'-sr^+ • ■ ■ +>''s^-'> 

O^Jf-P+2 OXif-p^2 O^^tr-P-^2 


I , 2»<pi , btps , , - bipp _- 

s:T" + + Xaj^ +...+ Xp— — 0 


when they are considered as linear homogeneous equations in the quantities 1, Xi, Xs, . . . Xp. 

Since a determinant (A7) exists for each value of M up to {N — P), a set of equations (A8) exists for each A/. It is 
seen that in these sets only the.first equation varies, since only tlie first equation depends upon M. Observation 
shows that all the {N — P) determinants of equation (A7) can be derived from the set of N equations 


^+Xi^‘+Xa 


— r • • ^ 


by successively writing the determinants of compatibility 
of the last P equations with each of the first (N—P) equa- 
tions in turn. It has thus been proven that if equations (A3) 
are true, the minimizing equations (A6), equivalent to 
equations (A7), must hold. 

It is seen, however, that equations (A3) are N equations 
in (JV+P) unkno%vns consisting of N »’s and P X’s. The 
remaining necessary P equations come from the original 
equations of constraint (A2). Hence, the simultaneous 
equations (A2) and (A3) 


+Xp|g=0 (K=1,2,3, 


<Pj{Zu Xi, Xg, . . . Xn)=0 (J=1, 2, 3, . . . F) (A2) 


are necessary equations for the minimization of J{xi, Xq, Xt, 
ijv), which was to be proved. 
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